Scale invariance of entanglement dynamics in Grover's quantum search algorithm 
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We calculate the amount of entanglement of the multiqubit quantum states employed in the 
Grover algorithm, by following its dynamics at each step of the computation. We show that genuine 
multipartite entanglement is always present. Remarkably, the dynamics of any type of entanglement 
as well as of genuine multipartite entanglement is independent of the number n of qubits for large 
n, thus exhibiting a scale invariance property. We compare this result with the entanglement 
dynamics induced by a fixed-point quantum search algorithm. We also investigate criteria for 
efficient simulatability in the context of Grover's algorithm. 
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I. INTRODUCTION 



Although it is well-known that entanglement repre- 
sents an essential ingredient in quantum communication 
and information, its role in the speed-up of quantum com- 
putational processes is not yet fully understood and still 
represents a debated question [TJH]. In particular, it is 
of great interest to investigate the role of multipartite 
entanglement in quantum algorithms. In Shor's algo- 
rithm multipartite entanglement was proved to be neces- 
sary to achieve exponential computational speed-up with 
quantum resources £Q. Moreover, more recently it was 
shown that multipartite entangled states are employed in 
the Deutsch-Jozsa algorithm and in the first step of the 
Grover algorithm [5j . In our work we investigate in detail 
the entanglement properties in the Grover algorithm [5], 
namely we study the behaviour of entanglement of the 
states of n qubits along the whole computational pro- 
cess, and disclose in particular a noteworthy scale invari- 
ance property of its dynamics in terms of the geometric 
measure of entanglement (GME). Previous works on the 
entanglement dynamics in Grover's algorithm considered 
other entanglement measures and focused only on bipar- 
tite entanglement (see, for example, 011]). We also study 
the entanglement dynamics in the fixed-point 7r/3 quan- 
tum search algorithm ]§\ [TU] and show that it turns out 
to be qualitatively similar to the Grover case. 

This paper is organised as follows. In Sect. [IT] we con- 
sider the Grover quantum search algorithm and study its 
entanglement dynamics for any number of qubits in the 
cases of one and two solutions to the search problem. In 
Sect. |III| we compare these results with the entangle- 
ment dynamics in fixed-point search algorithms. Finally, 
in Sect. II VI we summarise the main results and comment 
on their possible relations to classical efficient simulata- 
bility of Grover's algorithm. 



II. ENTANGLEMENT DYNAMICS IN 
GROVER'S ALGORITHM 

Let us remind the reader that the Grover search al- 
gorithm |B] employs pure states of n qubits which are 
initially prepared in an equally weighted superposition 
of all computational basis states \ipo) = ^J^q 1 \x), 
which can be more conveniently written as 
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where N = 2™ and M is the number of searched items 
(in the following also referred to as "solutions" of the 
search problem). Here, = — j=? Y] x \x s ) represents 
the superposition of all the states \x, 
(i.e. searched items), and \X ) = , , , , > , ,r 



that are solutions 

2rr I*") ^e- 



y/N-M 

notes the superposition of all the states \x n ) that are not 
searched for. The global state after k iterations of the 
Grover operation G has the form [TTJ [T^] 



l^fe) = G fe |Vo) = cos9 k \X ) +sin0 fe |X 1 ), (2) 



with 9 k = (k + 1/2)9 and 9 = 2^/M/N in the limit M < 
N. The unitary operation G is usually decomposed in 
two basic blocks, G — IU, where U represents the oracle 
call, i.e. U — 1 —2\Xi) (Xi |, and 1 is the inversion about 
the mean operation, namely X — —(1 — 2\i\jq) (ipo |). The 
operation G is repeated until the state \ipk) overlaps as 
much as possible with |-Xi), namely for k opt = CI[(n/9 — 
l)/2], where CI[x] denotes the closest integer to x. In 
the limit M < C N, t he optimal number of iterations is 
k op t = CI[jy/N/M — |], i.e. it is proportional to the 
square root of N. In the following we will consider the 
condition M <C N to be always fulfilled. 

We will now study the entanglement properties of the 
states ([2]) as functions of the number of iterations k and 
the number of qubits n for a fixed number of solutions. 
We will quantify the amount of entanglement by the 
GME [13 , which for a pure n-partite state reads 
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where S q is the set of g-separable states, namely states 
that are separable for q partitions of the n-qubit system. 
The GME represents a suitable entanglement measure 
when multi-partite systems are taken into account. No- 
tice that E n quantifies the amount of entanglement of 
any kind contained in the global system, i.e. it is non- 
vanishing even for states showing entanglement just be- 
tween two subsystems, while E 2 quantifies genuine mul- 
tipartite entanglement [14] . 



A. Single solution 

Let us first consider the case of a single solution to 
the search problem (M = 1). W.l.o.g., as will be proved 
later, we consider the state \X\) representing the solution 
to be invariant under any permutation of the n qubits 
(e.g. 1 1 1 1 ... 1) ) . Therefore, the state \ipk,M=i) at s * e P ^ 
of the algorithm is also permutation invariant for all fc's. 
Let us first compute E n for this set of states for vary- 
ing k. Due to this symmetry property, the search for the 
maximum in Eq. (T3J) can be restricted to symmetric sepa- 
rable states \d>)® n so that the maximisation involves 
only the two parameters a € [0, 7r] and f3 € [0,27r] that 
define the single qubit state \(j>) = cos f |0) + e I/3 sin §|1). 
Furthermore, since 9k <E [0,7r/2] the coefficients of \ipk) 
are all positive and the optimal value of the phase factor 
can be fixed to j3 = 0. 

The GME E n for a single solution then takes the form 
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(4) 



The optimal value of a can then be found by setting 
t = tan ^ and calculating the derivative of the overlap 
explicitly, which reduces to finding the root of a polyno- 
mial in t. 

In Fig. [Tlwe report the behaviour of E n {\i)jk,M=\)) for 
n = 12: The entanglement increases in the first half of 
iterations, achieves the maximal value of about 1/2, and 
then decreases to zero as soon as the optimal number 
of iterations is reached. This behaviour is qualitatively 
similar to the ones shown in 7, 8 , where the dynamics of 
both the two-qubit concurrence and the Von Neumann 
entropy of the half-qubit reduced state was studied. 

In order to quantify only genuine multipartite entan- 
glement we will now calculate E 2 . The expression of 
^2(1*0)) can be rewritten as [H] 



E2(\ip)) = 1 — max max /i 
p 1-1 



(5) 



where the /i's are the Schmidt coefficients of \tp) with re- 
spect to a fixed bipartition P : Q, and maxp denotes the 
maximisation over all possible bipartitions. Notice that, 
since the regarded state \ipk,M=i) is permutation invari- 
ant, we need to check only \n/2\ bipartitions, where \_x\ 
is the largest integer smaller or equal to x. In order to 
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Figure 1. (Color online) Evolution of entanglement as a func- 
tion of the number of steps k, for a single searched item, with 
n = 12 qubits and k opt = 49. E„(\ipk,M=i)) is depicted by 
blue dots, while E2(\ ipk,M=i}) by purple squares. The yellow 
dots represent the success probability. 



find the maximal Schmidt coefficient of \ipk.M=i) among 
all possible bipartitions we fix a generic bipartite split- 
ting P : Q, where P is composed of m qubits and Q of 
the remaining n — m, and compute the eigenvalues of the 
reduced density operator p P = TrQ[|-0 M/=1 )(V> fc ,M=i |], 
given by the following 2 m x 2 m matrix 



PP 



fa 



a 



a b 
b c 



(6) 



where a = 2 n ~ m A 2 , b = a - A(A - B), and c = a - 
A 2 + B 2 , with A = cos0 fe /V2" - 1 and B = sin6» fc . The 
maximal eigenvalue of the above matrix is given by 



.l)(2 n - m -I) A 2 {A- B) % Y . (7) 



The above expression shows that the bipartition that 
leads to the maximum eigenvalue corresponds to m = 1 
for all values of k. According to Eq. ([5|, the multipartite 
GME E 2 takes the explicit form 



^2(|^fe,M=l» = 



(8) 
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, cos 9 k 

' y/2 n - 1 
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This result shows that genuine multipartite entanglement 
has a qualitative similar behaviour as £ , n (|V'fe,M=i)) (see 
Fig. IJ, even if it achieves a maximum of about 0.14 
and the curve is derivable in that point. Notice also that 
-E'2(|V ; fe,M=i)) is symmetric with respect to k op t/2. 

We will now show that the entanglement dynamics in 
the Grover algorithm, namely the behaviour of E n and 
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E 2 during the operation of the algorithm, does not de- 
pend on the number of qubits n, thus exhibiting the prop- 
erty of scale invariance. For 2" S> 1 the two entanglement 
measures take the simple forms 



E n {\i>k,M=i))- 
E 2 {\^k.M=i))- 



sin 2 6> fe for 6> fe < tt/4, 
cos 2 9k for k > 7r/4, 



1 



sin^ 20 k 



(9) 



Both expressions depend only on 6k ~ ^k/k opt , namely 
on k/k pt, and not on k and n separately. Therefore, the 
entanglement dynamics of the Grover algorithm is scale 
invariant in the sense that it only depends on the number 
of steps taken, relative to the total number, but not on 
the length of the list. 

We want to point out that all the results presented 
so far, even if they were explicitly derived for permuta- 
tion invariant states, hold for any instance of the Grover 
algorithm with one searched item, i.e. M = 1. The 
number of possible single searched items in the Grover 
algorithm is 2™, which corresponds to the number of dis- 
tinct states All of these states can be achieved 
from a symmetric one by applying tensor products of a x 
Pauli operators and identity operators 1 (e.g. |001...1) = 
c x i ® a x2 <E> I3 Since these operations are lo- 
cal, they do not change the entanglement content of the 
resulting state. 



B. Two solutions 

Let us now consider the case of two searched items (i.e. 
M = 2). As an illustrative example we will consider the 
case in which both 1 00 . . . 0) and |11 . . . 1) arc solutions of 
the search problem, thus the state \Xi) is a GHZ state 
composed of n qubits, and the state at each step of the 
computation is permutation invariant. The measure of 
any entanglement E n is given by 



E n (\lpk,M=2)) = 1 - 



cos 9 k 



( 



cos™ I + sin" | 



V2™ - 2 
sin 9k 



a . a>i 
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We maximised this quantity numerically; in Fig. [2] we 
show the behavior for n — 13. Notice that after k opt iter- 
ations, the measure E n (\ipk M=2)) is no longer zero but 
equal to 1/2. That is because the final state is no longer 
fully separable but instead it is close to the GHZ state. 
In this case the maximal value that the entanglement 
reaches during the computation is about 2/3, i.e. higher 
than the case M = 1. Furthermore, this maximum is 
no longer reached at half of the optimal number of steps 
k pt, but in a later step, i.e. k/k op t — 0.61. 

Regarding genuine multipartite entanglement, E 2 with 
two symmetric solutions can be computed by following an 
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Figure 2. (Color online) The GME as a function of the number 
of steps k, for two symmetric solutions of the search problem. 
The number of qubits is n = 13, and k opt = 49. E n (\if)k,M=2}) 
is given by blue dots, E2(\tpk,M=2)) by purple squares. The 
yellow dots represent the success probability. 



analogous procedure as for a single solution. The reduced 
density matrix for the general bipartite splitting P : Q, 
where m qubits are in P and n — m in Q, is now given 
by 



/c b 
b a 



Pp 



b a 
\d b 



b d\ 

a b 



a b 
b cj 



(10) 



where a,b,c and A, B are given below Eq. (6|, and 
d = a — 2A(A — B). It turns out that again the max- 
imum eigenvalue corresponds to the bipartite splitting 
with m — 1, and E 2 (\'4 ! k,M=^)) can be expressed analyt- 
ically as 



E 2 (\ipk,M=z)) = 1- 



2 n - 4 
2 n - 2 



cos V k - 



cos Ok sin Ok 



V2" - 2 ^2 

(11) 

This result is shown in Fig. [2] Notice that mul- 
tipartite entanglement has a different behaviour from 
En(\ipk,M= 2 ))- It is a monotonically increasing function 
that approaches the maximum value of 1/2 when the 
computation stops. 

In the asymptotic limit 2" ^> 1 the GME can be ex- 
pressed as 



E n (\lpk,M=2))- 



E 2 {\^k,M=-2)) r - 



sin 2 9k 

l+cos 2 9 k 



sin Ok ■ 



for 9k < arccos 1 / \/3, 
for 9k > arccos 1 / v3 



(12) 



As a consequence, both quantities exhibit the same scale 
invariance behaviour as discussed above for the case with 
one searched item, i.e. M — 1. We point out that the 
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Figure 3. (Color online) Evolution of entanglement in the 
7r/3 search as a function of the recursive step m, for a single 
searched item, with n = 12 qubits. Here and below in Fig. [4] 
blue dots, purple squares and yellow dots denote E n , E2 and 
the success probability, respectively. 

above results can be generalized to those search problems 
in which the two solutions are different in all digits, but 
not to all search problems with M = 2. 

III. ENTANGLEMENT DYNAMICS IN THE 
FIXED-POINT tt/3 QUANTUM SEARCH 
ALGORITHM 

In the previous section we have shown that the amount 
of entanglement is non- vanishing during the Grover algo- 
rithm and that its behaviour is scale invariant for a single 
solution to the search problem and in some instances of 
two solutions. We will now show that a similar entan- 
glement dynamics can be found in the fixed-point tt/3 
quantum search. This kind of quantum search algorithm 
was first introduced in [9] to overcome the fact that the 
Grover algorithm might lead to a high error probability 
if the number of solutions M is unknown, since it re- 
quires to stop at a precise iteration k opt , which depends 
on M. In contrast the tt/3 quantum search always con- 
verges to the solutions, and thus it can be regarded as 
a fixed-point algorithm, even if it is never as fast as the 
standard Grover algorithm. 

A possible way to realise such a fixed-point search is 
to slightly modify the operations U and X in order to 
produce a tt/3 shift instead of a tt shift [9], i.e. 

U^U f =1-(1- e i ^)\X 1 )(X 1 1, (13) 

I— >I S =-(l-(l-e i *)|^)<Vto|). 

Then, the sequence of gates to be applied is defined by 
the following recursive formula 

= A m l*A^JJ % A m , (14) 
A = l. 
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Figure 4. (Color online) Evolution of entanglement in the n/3 
search as a function of the recursive step m, for two symmetric 
solutions of the search problem, with n = 13 qubits. 

We now compute both E n and E 2 for the employed 
states at each recursive step m of the evolution. The re- 
sults were obtained numerically and are shown in Figs. 
H] and 2] for both one and two solutions. Notice that the 
entanglement behaviour of both E n and E2 is similar to 
the dynamics of the standard Grover algorithm. These 
results indicate that entanglement plays the same crucial 
role in both algorithms, even if a scale invariance prop- 
erty cannot be proved in the fixed-point algorithm case. 

IV. CONCLUSIONS 

In summary, we have studied the evolution of entan- 
glement in Grover's algorithm (for a small number of 
searched items), quantifying it via the GME. In particu- 
lar, we give an explicit formula for the amount of genuine 
multipartite entanglement, which is proven to be always 
non-zero throughout the computation. Interestingly, the 
dynamics of entanglement shows the behaviour of scale 
invariance, i.e. counter-intuitively the amount of entan- 
glement employed in the algorithm does not depend on 
the length of the searched list, but only on the number 
of steps taken, relative to the optimal number of steps. 
Since scale invariance is an important phenomenon in 
several areas of physics and mathematics, our results may 
open new avenues in the understanding of scale invari- 
ance properties of entanglement in other contexts, such 
as for example in many-body systems and phase tran- 
sitions. We have also compared the Grover search en- 
tanglement dynamics with the one of a different kind of 
search algorithm, i.e. the tt/3 quantum search, and we 
have showed that they exhibit a similar behavior. 

As a final comment, we may wonder whether the 
presence of true multipartite entanglement means that 
Grover's algorithm cannot be simulated efficiently by 
classical means. By efficient classical simulation of 
Grover's algorithm we mean that, given a database of 
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n qubits, i.e. 2™ items, it is possible to classically sim- 
ulate Grover's algorithm with a total cost that scales as 
V / 2™poly(n). We will now show that well-known crite- 
ria which guarantee efficient simulatability do not apply 
for Grover's algorithm. According to the Gottesman- 
Knill theorem [2] [12] , if a quantum computation starts 
in a computational basis state and involves only stabilizer 
gates (i.e. Hadamard, CNOT, phase gates and measure- 
ment of operators in the Pauli group), then it can be 
efficiently simulated on a classical computer. However, it 
can be easily shown that T transforms an element belong- 
ing to the Pauli group, i.e. a z (g) D™ -1 ), to an operator 
that no longer belongs to the Pauli group, and there- 
fore it cannot be implemented by stabilizer gates. Let 
us also consider the simulatability criterion introduced 
in [3J, based on the maximal Schmidt rank \ of IV 1 ) over 
all possible bipartitions. According to [3J, if \ does not 
exceed poly(n) in a computation that consists of poly(n) 
elementary gates (i.e. one- and two-qubit gates), then 
the computation can be classically simulated efficiently. 
We notice that for states of the form ([2]), x is upper 
bounded by M + 1. However, although there exists a 



decomposition of the Grover operation G into poly(n) 
elementary gates O [17], the state after the action of 
each two-qubit gate does not have a simple symmetric 
form and we no longer can keep track of the maximal 
Schmidt rank. Therefore we cannot conclude efficient 
simulatability. 

The above results show that, although the Grover op- 
eration cannot be implemented by stabilizer gates and 
therefore the Knill-Gottesman theorem cannot be ap- 
plied, the employed states at each Grover iteration are 
only slightly entangled according to the criterion sug- 
gested in [3]. These insights are nevertheless not suf- 
ficient to answer the question of simulatability of the 
Grover algorithm, which at present remains open. 
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